Effect of exercise on fluoride metabolism in adult humans: a pilot study

An understanding of all aspects of fluoride metabolism is critical to identify its biological effects and avoid fluoride toxicity in humans. Fluoride metabolism and subsequently its body retention may be affected by physiological responses to acute exercise. This pilot study investigated the effect of exercise on plasma fluoride concentration, urinary fluoride excretion and fluoride renal clearance following no exercise and three exercise intensity conditions in nine healthy adults after taking a 1-mg Fluoride tablet. After no, light, moderate and vigorous exercise, respectively, the mean (SD) baseline-adjusted i) plasma fluoride concentration was 9.6(6.3), 11.4(6.3), 15.6(7.7) and 14.9(10.0) ng/ml; ii) rate of urinary fluoride excretion over 0–8 h was 46(15), 44(22), 34(17) and 36(17) μg/h; and iii) rate of fluoride renal clearance was 26.5(9.0), 27.2(30.4), 13.1(20.4) and 18.3(34.9) ml/min. The observed trend of a rise in plasma fluoride concentration and decline in rate of fluoride renal clearance with increasing exercise intensity needs to be investigated in a larger trial. This study, which provides the first data on the effect of exercise with different intensities on fluoride metabolism in humans, informs sample size planning for any subsequent definitive trial, by providing a robust estimate of the variability of the effect

Capturing the essence of folding and functions of biomolecules using Coarse-Grained Models

The distances over which biological molecules and their complexes can function range from a few nanometres, in the case of folded structures, to millimetres, for example during chromosome organization. Describing phenomena that cover such diverse length, and also time scales, requires models that capture the underlying physics for the particular length scale of interest. Theoretical ideas, in particular, concepts from polymer physics, have guided the development of coarse-grained models to study folding of DNA, RNA, and proteins. More recently, such models and their variants have been applied to the functions of biological nanomachines. Simulations using coarse-grained models are now poised to address a wide range of problems in biology.Comment: 37 pages, 8 figure

The Origin of Minus-end Directionality and Mechanochemistry of Ncd Motors

Adaptation of molecular structure to the ligand chemistry and interaction with the cytoskeletal filament are key to understanding the mechanochemistry of molecular motors. Despite the striking structural similarity with kinesin-1, which moves towards plus-end, Ncd motors exhibit minus-end directionality on microtubules (MTs). Here, by employing a structure-based model of protein folding, we show that a simple repositioning of the neck-helix makes the dynamics of Ncd non-processive and minus-end directed as opposed to kinesin-1. Our computational model shows that Ncd in solution can have both symmetric and asymmetric conformations with disparate ADP binding affinity, also revealing that there is a strong correlation between distortion of motor head and decrease in ADP binding affinity in the asymmetric state. The nucleotide (NT) free-ADP (?-ADP) state bound to MTs favors the symmetric conformation whose coiled-coil stalk points to the plus-end. Upon ATP binding, an enhanced flexibility near the head-neck junction region, which we have identified as the important structural element for directional motility, leads to reorienting the coiled-coil stalk towards the minus-end by stabilizing the asymmetric conformation. The minus-end directionality of the Ncd motor is a remarkable example that demonstrates how motor proteins in the kinesin superfamily diversify their functions by simply rearranging the structural elements peripheral to the catalytic motor head domain

A Computational Investigation on the Connection between Dynamics Properties of Ribosomal Proteins and Ribosome Assembly

Assembly of the ribosome from its protein and RNA constituents has been studied extensively over the past 50 years, and experimental evidence suggests that prokaryotic ribosomal proteins undergo conformational changes during assembly. However, to date, no studies have attempted to elucidate these conformational changes. The present work utilizes computational methods to analyze protein dynamics and to investigate the linkage between dynamics and binding of these proteins during the assembly of the ribosome. Ribosomal proteins are known to be positively charged and we find the percentage of positive residues in r-proteins to be about twice that of the average protein: Lys+Arg is 18.7% for E. coli and 21.2% for T. thermophilus. Also, positive residues constitute a large proportion of RNA contacting residues: 39% for E. coli and 46% for T. thermophilus. This affirms the known importance of charge-charge interactions in the assembly of the ribosome. We studied the dynamics of three primary proteins from E. coli and T. thermophilus 30S subunits that bind early in the assembly (S15, S17, and S20) with atomic molecular dynamic simulations, followed by a study of all r-proteins using elastic network models. Molecular dynamics simulations show that solvent-exposed proteins (S15 and S17) tend to adopt more stable solution conformations than an RNA-embedded protein (S20). We also find protein residues that contact the 16S rRNA are generally more mobile in comparison with the other residues. This is because there is a larger proportion of contacting residues located in flexible loop regions. By the use of elastic network models, which are computationally more efficient, we show that this trend holds for most of the 30S r-proteins

Exploring the Conformational Transitions of Biomolecular Systems Using a Simple Two-State Anisotropic Network Model

Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called ANMPathway, for constructing a physically reasonable pathway between two endpoints of a conformational transition. We adopt a coarse-grained representation of the protein and construct a two-state potential by combining two elastic network models (ENMs) representative of the experimental structures resolved for the endpoints. The two-state potential has a cusp hypersurface in the configuration space where the energies from both the ENMs are equal. We first search for the minimum energy structure on the cusp hypersurface and then treat it as the transition state. The continuous pathway is subsequently constructed by following the steepest descent energy minimization trajectories starting from the transition state on each side of the cusp hypersurface. Application to several systems of broad biological interest such as adenylate kinase, ATP-driven calcium pump SERCA, leucine transporter and glutamate transporter shows that ANMPathway yields results in good agreement with those from other similar methods and with data obtained from all-atom molecular dynamics simulations, in support of the utility of this simple and efficient approach. Notably the method provides experimentally testable predictions, including the formation of non-native contacts during the transition which we were able to detect in two of the systems we studied. An open-access web server has been created to deliver ANMPathway results. © 2014 Das et al

A-Site Residues Move Independently from P-Site Residues in all-Atom Molecular Dynamics Simulations of the 70S Bacterial Ribosome

The ribosome is a large macromolecular machine, and correlated motion between residues is necessary for coordinating function across multiple protein and RNA chains. We ran two all-atom, explicit solvent molecular dynamics simulations of the bacterial ribosome and calculated correlated motion between residue pairs by using mutual information. Because of the short timescales of our simulation (ns), we expect that dynamics are largely local fluctuations around the crystal structure. We hypothesize that residues that show coupled dynamics are functionally related, even on longer timescales. We validate our model by showing that crystallographic B-factors correlate well with the entropy calculated as part of our mutual information calculations. We reveal that A-site residues move relatively independently from P-site residues, effectively insulating A-site functions from P-site functions during translation

Shrub Invasion Decreases Diversity and Alters Community Stability in Northern Chihuahuan Desert Plant Communities

BACKGROUND:Global climate change is rapidly altering species range distributions and interactions within communities. As ranges expand, invading species change interactions in communities which may reduce stability, a mechanism known to affect biodiversity. In aridland ecosystems worldwide, the range of native shrubs is expanding as they invade and replace native grassland vegetation with significant consequences for biodiversity and ecosystem functioning. METHODOLOGY:We used two long-term data sets to determine the effects of shrub encroachment by Larrea tridentata on subdominant community composition and stability in formerly native perennial grassland dominated by Bouteloua eriopoda in New Mexico, USA. PRINCIPAL FINDINGS:Our results indicated that Larrea invasion decreased species richness during the last 100 years. We also found that over shorter temporal scales species-poor subdominant communities in areas invaded by Larrea were less stable (more variable in time) compared to species rich communities in grass-dominated vegetation. Compositional stability increased as cover of Bouteloua increased and decreased as cover of Larrea increased. SIGNIFICANCE:Changes in community stability due to altered interspecific interactions may be one mechanism by which biodiversity declines in grasslands following shrub invasion. As global warming increases, shrub encroachment into native grasslands worldwide will continue to alter species interactions and community stability both of which may lead to a decline in biodiversity

Multiple Routes and Milestones in the Folding of HIV–1 Protease Monomer

Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV–1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV–1 protease

Minimum Free Energy Path of Ligand-Induced Transition in Adenylate Kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme

Interconversion of Functional Motions between Mesophilic and Thermophilic Adenylate Kinases

Dynamic properties are functionally important in many proteins, including the enzyme adenylate kinase (AK), for which the open/closed transition limits the rate of catalytic turnover. Here, we compare our previously published coarse-grained (double-well Gō) simulation of mesophilic AK from E. coli (AKmeso) to simulations of thermophilic AK from Aquifex aeolicus (AKthermo). In AKthermo, as with AKmeso, the LID domain prefers to close before the NMP domain in the presence of ligand, but LID rigid-body flexibility in the open (O) ensemble decreases significantly. Backbone foldedness in O and/or transition state (TS) ensembles increases significantly relative to AKmeso in some interdomain backbone hinges and within LID. In contact space, the TS of AKthermo has fewer contacts at the CORE-LID interface but a stronger contact network surrounding the CORE-NMP interface than the TS of AKmeso. A “heated” simulation of AKthermo at 375K slightly increases LID rigid-body flexibility in accordance with the “corresponding states” hypothesis. Furthermore, while computational mutation of 7 prolines in AKthermo to their AKmeso counterparts produces similar small perturbations, mutation of these sites, especially positions 8 and 155, to glycine is required to achieve LID rigid-body flexibility and hinge flexibilities comparable to AKmeso. Mutating the 7 sites to proline in AKmeso reduces some hinges' flexibilities, especially hinge 2, but does not reduce LID rigid-body flexibility, suggesting that these two types of motion are decoupled in AKmeso. In conclusion, our results suggest that hinge flexibility and global functional motions alike are correlated with but not exclusively determined by the hinge residues. This mutational framework can inform the rational design of functionally important flexibility and allostery in other proteins toward engineering novel biochemical pathways